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Abstract 

Cartesian products of graphs and hypergraphs have been studied since the 1960s. For (un)directed hypergraphs, 
unique prime factor decomposition (PFD) results with respect to the Cartesian product are known. However, there is 
still a lack of algorithms, that compute the PFD of directed hypergraphs with respect to the Cartesian product. 

In this contribution, we focus on the algorithmic aspects for determining the Cartesian prime factors of a finite, 
connected, directed hypergraph and present a first polynomial time algorithm to compute its PFD. In particular, the 
algorithm has time complexity 0{\E\\V\r~) for hypergraphs H - (V, £), where the rank r is the maximum number 
of vertices contained in an hyperedge of H. If r is bounded, then this algorithm performs even in OdZsl log^(|y|)) 
time. Thus, our method additionally improves also the time complexity of PFD-algorithms designed for undirected 
hypergraphs that have time complexity Ddfiliyir^A®), where A is the maximum number of hyperedges a vertex is 
contained in. 

Keywords: Directed Hypergraph, Cartesian Product, Prime Factor Decomposition, Factorization Algorithm, 
2-Section 


1. Introduction 


Products are a common way in mathematics of constructing larger objects from smaller building blocks. For 
graphs, hypergraphs, and related set systems several types of products have been investigated, see lll8lll4ll for recent 
overviews. 

In this contribution we will focus on the Cartesian product of directed hypergraphs that are the common gen¬ 
eralization of both directed graphs and (undirected) hypergraphs. In particular, we present a fast and conceptually 
very simple algorithm to find the decomposition of directed hypergraphs into prime hypergraphs (its so-called prime 
factors), where a (hyper)graph is called prime if it cannot be presented as the product of two nontrivial (hyper)graphs, 
that is, as the product of two (hyper)graphs with at least two vertices. 

Graphs and the Cartesian Product. A graph is a tuple G = {V, E) with non-empty set of vertices V and a set of 
edges E containing two-element subsets of V. If the edges are ordered pairs, then G is called directed and undirected, 
otherwise. The Cartesian graph product was introduced by Gert Sabidussi lE^I . As noted by Szamkolowicz 
also Shapiro introduced a notion of Cartesian products of graphs in 12711 . Sabidussi and independently V.G. Vizing 
S showed that connected undirected graphs have a representation as the Cartesian product of prime graphs that is 
unique up to the order and isomorphisms of the factors. The question whether one can find the prime factorization 
of connected undirected graphs in polynomial time was answered about two decades later by Feigenbaum et al. in 
who presented an (9(|y|^-^) time algorithm. From then on, a couple of factorization algorithms for undirected graphs 
have been developed |[l| [Ill [o , U ^ . The fastest one is due to Imrich and Peterin 121 ] and runs in linear-time 

om + i£i). 

For connected directed graphs, Feigenbaum showed that the Cartesian product satisfies the unique prime factor¬ 
ization property d. Additionally, she provided a polynomial-time algorithm to determine the prime factors which 
was improved to a linear time approach by Crespelle et al. (@]. 
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Hypergraphs and the Cartesian Product. Hypergraphs are a natural generalization of graphs in which “edges” may 
consist of more than two vertices. More precisely, a hypergraph is a tuple H - {V, E) with non-empty set of vertices 
V and a set of hyperedges E, where each e e £ is an ordered pair of non-empty sets of vertices e = it(e), h{e)). If 
t{e) = h{e) for all e e £ the hypergraph is called undirected and directed, otherwise. Products of hypergraphs have 


13, [3, in 0,111 ms mill. It was shown 


been investigated by several authors since the 1960s 

by Imrich il9l] that connected undirected hypergraphs have a unique prime factor decomposition (PFD) w.r.t. to the 
Cartesian product, up to isomorphism and the order of the factors. A first polynomial-time factorization algorithm for 
undirected hypergraphs was proposed by Bretto et al. H. 

Unique prime factorization properties for directed hypergraphs were derived by Ostermeier et al. 10]. However, 
up to our knowledge, no algorithm to determine the Cartesian prime factors of a connected directed hypergraph is 
established, so-far. 

Summary of the Results. In this contribution, we present an algorithm to compute the PFD of connected directed hy¬ 
pergraphs in 0{\V\\E\r^) time, where the rank r denotes the maximum number of vertices contained in the hyperedges. 
In addition, if we assume to have hypergraphs with bounded rank the algorithm runs in 0{\E\ log^dV])) time. Note, as 
directed hypergraphs are a natural generalization of undirected hypergraphs, our method generalizes and significantly 
improves the time-complexity of the method by Bretto et al. |@]. In fact, the algorithm of Bretto et al. has time com¬ 
plexity (9(|y||£|A®r®), where A is the maximum number of hyperedges a vertex is contained in. Assuming that given 
hypergraphs have bounded rank r and bounded maximum degree A this algorithm runs therefore in (9(|y||£|) time. 

We shortly outline our method. Given an arbitrary connected directed hypergraph H - (y £) we first compute 
its so-called 2-section [//] 2 , that is, roughly spoken the underlying undirected graph of H. This allows us to use the 
algorithm of Imrich and Peterin 112111 in order to compute the PFD of [H ]2 w.r.t. the Cartesian graph product. As we 
will show, this provides enough information to compute the Cartesian prime factors of the directed hypergraph H. In 
distinction from the method of Bretto et al. our algorithm is in a sense conceptually simpler, as (7 j we do not need the 
transformation of the hypergraph H into its so-called L2-section and back, where the L2-section is is an edge-labeled 
version of the 2-section [7/]2, and (2) the test which (collections) of the factors of the 2-section are prime factors 
of H follows a complete new idea based on increments of fixed vertex-coordinate positions, that allows an easy and 
efficient check to determine the PFD of H. 


2 . Preliminaries 

2.1. Basic Definitions 

A directed hypergraph H - (V,E) consists of a finite vertex set V{H) := V and a set of directed hyperedges or 
(hyper)arcs E{H) ;= £. Each arc e e £ is an ordered pair of non-empty sets of vertices e = (t(e),h{e)). The sets 
t(e) c y and h{e) c y are called the tail and head of e, respectively. The set of vertices, that are contained in an arc 
will be denoted by V(e) t(e)Uh(e). If t(e) - h(e) holds for all e e £, we identify e with V(e), and we call H - (V, E) 
an undirected hypergraph. An undirected hypergraph is an undirected graph if |y(e))| = 2 for all e e E. The elements 
of £ are called simply edges, if we consider an undirected graph. The hypergraph with jyj = 1 and £ = 0 is denoted 
by 7^1 and is called trivial. 

Throughout this contribution, we only consider hypergraphs without multiple hyperedges and thus, being £ a 
usual set, and without loops, that is, |y(e)| > 1 holds for all e e £. However, we allow to have hyperedges being 
properly contained in other ones, i.e., we might have arcs e, / 6 £ with t{e) c t{f) and h{e) c h{f). 

A partial hypergraph or sub-hypergraph H’ - (V', £') of a hypergraph H = (y £), denoted by 77' c 77, is a 
hypergraph such that V' C V and £' c £. The partial hypergraph 77' = (V',E') is induced (by V) if £' = {e 6 £ j 
y(e) c y'}. Induced hypergraphs will be denoted by (y'). 

A weak path P (joining the vertices VQ,Vk e y) in a hypergraph 77 = (K £) is a sequence P - 
(vo, ei, vi, 62 , ■ • ■, 6 t:, Vk) of distinct vertices and arcs of 77, such that vo € y(ei), Vk e Viek) and vj 6 Vief) n V(ej+\). A 
hypergraph 77 is said to be weakly connected or simply connected for short, if any two vertices of 77 can be joined by 
a weak path. A connected component of a hypergraph 77 is a connected sub-hypergraph H' c H that is maximal w.r.t. 
inclusion, i.e., there is no other connected sub-hypergraph 77" c 77 with 77' C 77". Usually, we identify connected 
components 77' = (V', £') of 77 simply by their vertex set V', since {V’) - 77'. 
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A homomorphism from Hi - (Vi,Ei) into H 2 - (^ 2 , £ 2 ) is a mapping (f> : Vi V 2 such that ^(e) is an arc in H 2 
whenever e is an arc in Hi with the property that (^(f(e)) = t((p{e)) and (p{h{e)) - h{(p(e)). A bijective homomorphism 
(f) whose inverse function is also a homomorphism is called an isomorphism. 

The rank of a hypergraph H - {V, E) is r(H) - maXggg \V(e)\. 

The 2-section {H '\2 of a (directed) hypergraph H - {V,E) is the undirected graph {V,E') with E' - 
{xy := {x,y} c y 1 2 c y, 3 e e £ : x,y c V{e)}. In other words, two vertices are linked by an edge in [//]2 if they 
belong to the same hyperarc in H. Thus, every arc e e £ of // is a complete graph in [//] 2 , i e-, all pairwise different 
vertices in V{e) are linked by an edge in [//] 2 - Complete graphs defined on a vertex set V will be denoted by K\v\. 

We will also deal with equivalence relations, for which the following notations are needed. For an equivalence 
relations R we write p = £ to indicate that g is an equivalence class of R. A relation Q is finer than a relation R while 
the relation R is coarser than Q if (e,/) e Q implies (e,/) e R, i.e, 2 c £. In other words, for each class p of £ there 
is a collection {x\x £ ^?) of Q-classes, whose union equals g. Equivalently, for all ^ = Q and i/^ = £ we have either 
(pCifforipnilf-d). 

Remark 1. If not stated differently, we assume that the hypergraphs considered in this contribution are connected. 


2.2. The Cartesian Product, (Pre-)Coordinates and (Pre-)Layers 

Let Hi and H 2 be two hypergraphs. The Cartesian product H - HiC\H 2 has vertex set V{H) - V{Hi) x V{H 2 ), 
that is the Cartesian set product of the vertex sets of the factors and the arc set 

E{H) ={({x! X t{f), (x) X h{f)) I X e V{Hi),f e £(// 2 )}U 
{(t(e) X {y} , h{e) x {y}) | e e E(Hi),y e V(H 2 )]. 

Thus, the tuple ({(x,,y,),/ e /}, {(x',yp ,7 e 7}) withx,, x' e V(Hi), yi,y'j e V(H 2 ), i e I,j e J is an arc in 
E(HiDH 2 ) if and only if either 


(i) ({x;, i 6 /}, {x'j, j e J]) is an arc in E{Hi) and y, = y' for all i e I,j e J or 

(ii) ({y,-, i e /}, {y', j e J]) is an arc in E{H 2 ) and x, = x' for all i e I,j e J. 


The Cartesian product is associative, commutative, and the trivial one-vertex hypergraph £1 without arcs serves 
as unit 118, 23] . Thus, for arbitrary finitely many factors {//, , i e /} the product □,£///, is well-defined, and each vertex 
X e y(n,e///,) is properly “coordinatized” by the vector (x,)fe/ whose entries are the vertices x; of the factors //,. 

A nontrivial hypergraph H is prime with respect to the Cartesian product if it cannot be represented as the Cartesian 
product of two nontrivial hypergraphs. A prime factor decomposition (PFD) of // is a representation as a Cartesian 
product H - □,£///, such that all factors i e /, are prime and Hi ^ Ki. Note, the number k of prime factors of 
H - (y £) is bounded by log{\V\), since every Cartesian product of k non-trivial hypergraphs has at least 2* vertices. 

Two important results concerning the Cartesian products of hypergraphs are given now. 


Lemma 2.1 The Cartesian product H - of directed hypergraphs is connected if and only if all of its 

factors Hi are connected. 


Theorem 2.2(H). Connected (directed) hypergraphs have a unique prime factor decomposition with respect to the 
Cartesian product. 


We will show, that the PFD of a hypergraph H can be obtained from the PFD of its 2-section [H] 2 - For this the 
following lemma is crucial. 

Lemma 2.3. is an arbitrary factorization of H it holds that [H )2 - '^iei[EIi\ 2 - 

Proof. Since the Cartesian product is commutative and associative it suffices to prove the statement for two factors. 
Assume that H - HiC\H 2 and every vertex x has coordinates (xi,X 2 ). Thus, there is an isomorphism f '■ V{H) —» 
V(HiOH 2 ) via X (xi, X 2 ). We show that f is also an isomorphism for the graphs {H '\2 and {Hi\ 2 ^{H 2 } 2 - 

The edge xy is contained in £([//] 2 ) if and only if there is an arc e e E{H) with x,y e V{e) - t{e) U h{e) if 
and only if (i) xi = yi and X 2 y 2 c V{f) e £( 1 / 2 ) or (ii) X 2 = y 2 and xiyi c V{f) e E(Hi) if and only if (i) 
xi = yi and X 2 y 2 € E([H 2 ] 2 ) or (ii) X 2 = y 2 and xiyi 6 E{{Hi\ 2 ) if and only if the edge (p(x)<p(y) is contained in 

£([//i]2n[//2]2)- □ 
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Now, given the PFD of H - □,£///,, we can infer that [H]2 - However, the factors [Hi\2 might 

not be prime w.r.t. the Cartesian graph product and hence, {H '\2 might have more prime factors. Since the PFD of 
{Hh — 'Oij^jGj is unique it follows that the 2-section of the prime factors Hi of // is a combination of the prime factors 
Gj of [H]2, that is, [//,]2 - J' c J for all i e I. 

Our algorithm will start with the PFD of [H ]2 = OjejGj w.r.t. the Cartesian product of undirected graphs and then 
tries to combine the respective prime factors of [//]2 to reconstruct the prime factors of H. In other words, we need 
to find suitable subsets J' c J so that [//,]2 - Ojej’Gj and H, is a prime factor of H. To this end, we will introduce 
(pre-)coordinates and (pre)-layers. 

Definition 2.4 ((Factorization) Coordinatization). Let H be isomorphic to some product □,£///„ where each factor Hj 
has vertex set {1,If. A factorization coordinatization or coordinatization/or short, is an isomorphism T from H 
to Thus, T assigns to a vertex v 6 V(H) a vector of coordinates {vfiei where 1 < V; < // is a vertex in V{Hi) 

Hence, a coordinatization gives in an explicit way the information of the underlying product structure of H. 
Hence, to find a factorization of H one can equivalently ask for a coordinatization of H, a fact that we will utilize 
in our algorithm. Note that the coordinatization w.r.t. a given product decomposition is unique up to relabeling the 
vertices in each factor //,. 

We will also need a notion which is similar to a coordinatization but is implied by a factorization of the 2-section 
[H ]2 rather than a decomposition of H. 

Definition 2.5 (Pre-Coordinatization). Let H be a given hypergraph and assume that [//]2 has a coordinatization 
T: y([//] 2 ) —> X,e/{1,..., If. Since VifH^f) — V(H) we infer that T is a bijective map on V(H) that assigns to each 
vertex v e V{H) a unique coordinate-vector (vfiei where 1 < v; < This map is called pre-coordinatization o///. 

For convenience, we will usually omit the function Y and identify every vertex v e V{H) with its (pre-)coordinate 
vector, i.e., we will write v = (v',■),■(=/ rather than T(v') = (v,■),£/. 

Definition 2.6 (Layers and Pre-Layers). Let H ^ □,£///, with given respective coordinatization v - (vfiei e V(H) 
and r c 1. The /'-layer through v (denoted by H'j,) with respect to this coordinatization is the sub-hypergraph induced 
by the vertices [u — (ufi^; \ i ( I' «,■ = v, }, i.e., we fix all coordinates except those contained in the set /'. Note, 

H'j, - 

Analogously, as layers are defined by means of a coordinatization define the pre-layers by means of a pre- 
coordinatization. 

For simplicity we write HJ instead ofHJ.^ and i-(pre-)layer rather than {i]-(pre-)layers. 

For later reference, we need the following observation and lemma. If // ^ Qie//// and e is an arc of H, then all 
vertices in V{e) are contained in the same Flayer H'.' for some i e I and v e V{H), i.e. they only differ in the Fth 
coordinate. The same is true for pre-layers. 

Lemma 2.7. Let H be a hypergraph and let T be a pre-coordinatization ofH. Then every arc e ofH contains vertices 
of exactly one pre-layer w.r.t T, that is, all vertices in V{e) only differ in the same i-th coordinate. 


Proof. Every hyperarc e forms a complete subgraph K\v(e)\ in [//] 2 - Moreover, complete subgraphs must be contained 
entirely in one of the /-layers of [//] 2 , as complete graphs are so-called S-prime graphs, see e.g. |lj,[iSll7ll2^. Hence, 
each hyperarc is contained in one /-pre-layer of H. □ 


Note, any isomorphism from [H ]2 to □,£/[//,]2 and thus, a pre-coordinatization of H, is a coordinatization of H if 
and only if H has a factorization □,£///,. Lemma 1231 immediately implies that every coordinatization of H is also a 
pre-coordinatization of H, while the converse is not true in general. On the other hand, we have the following result 
for so-called increments of coordinates. 


Definition 2.8 (Increments of Coordinates). Given a pre-coordinatization T; V{H) —> X,e/{1, ..., of H and a 
vertex v - (vi,..., v,■,..., v,t) 6 V{H) we define inc{v, i) (w.r.t. T) as the vertex with coordinates (vq,..., v',- -I- 1,..., Vk) 
where we set v',- -H 1 ;= 1 if Vi H- 1 > //. 

For an (ordered) set of vertices W <zV we define the (ordered) set inc(W, /)) = {inc(w, /) | w e W). 

Finally, we denote for an arc e — (t{e), h(e)) its increment (inc(t(e), /), inc{h{e), /)) by inc{e, /). 
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Lemma 2.9. Let H — {V, E) be a hypergraph and T: V —> X/(=/{l,..., /,} be a pre-coordinatization of H. If for each 
arc e e E (where the vertices of e differ only in the j-th coordinate), there is an arc inc{e, i) € E for all i + j, then T 
is a coordinatization of H. 

Proof By Lemma IZTl all vertices within one arc e 6 £ differ in precisely one coordinate. Let e e £ be an arbitrary 
hyperarc and assume the vertices differ in the j-th coordinate. 

Let Elj be the set of j-pre-layers contained in H. Let i + j be an arbitrary index i e I. Assume that for each 
hyperedge e contained in some j-pre-layer H{1) all “incremental copies” inc(e, i) are also contained in £, then there 
is a homomorphism from H{Y) - {V{H{\))) to H{2) - {{inc{V(H{l)), i))), where H{2) corresponds to some other 
j-pre-layer. Assume that for all such “consecutive” j-pre-layer there is a homomorphism from H{1) to H(l + 1), 
1 </</,- 1. By construction, after f - 1 incremental steps we arrive at the /,-th j-pre-layer //(/,) and hence, 
H(l) = ({inc(V(H{li)), /))). If there is an homomorphism from //(/,) to H{1), then there is trivially an isomorphism 
between all such j-pre-layers H{1),H{2),.. .,H{U) 6 “K/. Thus, if for all arcs e € £, where the vertices of e differ 
precisely in this j-th coordinate, there is a hyperarc inc(e, i) e £ for all i + j, then there isomorphism between all 
j-pre-layers contained in Eij for this fixed j 6 /. 

If this is true for all arcs e e £, and thus, for all /-pre-layers with i 6 /, then all such /-pre-layers are isomorphic 
for each i 6 /. 

In particular, we can define for vertices v, w and an index i e I the map gf^: H'j —> H'J which maps every vertex 
in H'.' to the unique vertex in H" with the same /-coordinate. By the preceding arguments, for each / € J the map g™ 
is an isomorphism between the /-pre-layers in H for all v, w e V(H). 

Finally, assume for contradiction that T is a not a coordinatization and hence, T is not an isomorphism from H 
to any product □,£///, . Hence, there must be some / 6 I such that not all /-layers are isomorphic by means of g™, a 
contradiction. □ 

Lemma lZTl allows defining an equivalence relation Rr on the hyperedge set E{H) for a given pre-coordinatization 
T: V(H) —> X,£/{1,..., /,} of H, as follows: (e,/) e Rr if e e H'i and / e Hf for some / e I and v, w e V(H). In other 
words, e and / are in relation Rr if they are both contained in the /-pre-layers for the same fixed / E I. Note, in case 
that T is a coordinatization the relation Rr is also known as product relation, that is, each equivalence of Rr contains 
the hyperedges of all copies of some (not necessarily prime) factor of H. In order to avoid confusion, we sometimes 
write that Rr{H) to indicate that Rr is defined on the edge set of H. 

Given two pre-coordinatizations Ti and T 2 , we say that Ti is finer than T 2 , while T 2 coarser than Ti if Rr^ is 
finer than Rr^- We can immediately infer the next result. 

Lemma 2.10. Let H be a hypergraph and let Ti and T 2 be pre-coordinatizations of H. Then Ti is finer than T 2 if 
and only if the factorization of[H ]2 corresponding to T 2 can be obtained from the factorization □,e/G; corresponding 
to Ti by combining some of these factors, i.e., I can be partitioned into I\,... ,Ii so that the pre-coordinatization T 2 
is an isomporphism between [H ]2 and where G'j ^ □fe/^.G;. 

Proof By definition T 1 is finer than T 2 if and only if Rr^ is finer than £x 2 • By definition of the Cartesian product and 
by Theorem l2.2l this is the case if and only if every factor in the factorization corresponding to T 2 is a combination of 
the factors in the coordinatization Ti. □ 

Remark 2. By Lemma 12.31 every coordinatization of H is also a pre-coordinatization of H. This implies together 
with Lemma \2. 1 OU hat every pre-coordinatization ofH can be achieved by starting with the pre-coordinatization w.r.t. 
the prime factorization of[H ]2 and then combining the corresponding pre-layers of H to obtain the layers w.r.t. the 
prime factorization of H. In other words, one needs to find a partition l+J^ey Ij of the index set I and then combine all 
pre-layers corresponding to indices in the same part into one. 

As we shall see later, in our algorithm we will only check increments w.r.t. the pre-coordinatization T' coming 
from the PFD of [//] 2 . However, we have to prove that this is indeed sufficient (Theorem 13. lb . In order to apply 
Lemma |2^ to validate whether we end up with a coordinatization T of // we would need to check increments with 
respect to this coarser (pre-)coordinatization T. Now one might hope that increments with respect to the coarser 
pre-coordinatization are automatically increments with respect to the finer pre-coordinatization T' or that at least the 
coarser pre-coordinates can be chosen in a suitable way. However, this is not the case as the following example shows. 
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Assume that at some point we need to combine i- and y'-pre-layers of sizes 3 and 4 respectively. The resulting 
A:-pre-layers with respect to the new coordinatization will each contain 12 vertices labeled 1,..., 12. We now claim 
that no matter how we assign the new labels, there is always at least one increment inc(-, k) which is not an increment 
inc(-, i) or inc(-, j). Assume for a contradiction that all increments inc(-, k) were either of the form inc(-, i) or inc(-, j). 
By applying inc(-, k) recursively 12 times to a vertex, we end up at the same vertex again. This means, that we have 
applied inc(-, i) a number of times which must be divisible by 3 and inc(-, j) a number of times which must be divisible 
by 4. However, no suitable multiples of 3 and 4 add up to 12. 

The latter example shows that the single check of increments with respect to the PFD of [//]2 is not sufficient to 
invoke Lemma |2^ to conclude that some coarser pre-coordinatization is indeed a coordinatization. For this purpose, 
we need the following additional lemma. 

Lemma 2.11. Let H — {V,E) be a hypergraph, letTi, T 2 be pre-coordinatizations ofH such that Ti is finer than T 2 . 
Let inc\ and inc 2 the respective increment maps. Assume that for each arc e e E (where the vertices of e differ only in 
the j-th coordinate w.r.t. T 2 ), there is an arc inc\{e, i) e E for all i i Ij (where Ij is defined as in Lemma [2.10\l . Then 
there is an arc inc 2 {e, k) for all k ^ j and hence T 2 is a coordinatization ofH. 

Proof. The coordinates of the vertices in inc 2 (e, k) w.r.t. Ti can be obtained from those of vertices in e by only 
changing coordinates outside Ij. This can be achieved by successive applications of inci(-,0 for i i Ij- Since we 
started at an edge (namely e) and each of those applications takes edges to edges we also end at an edge inc 2 (e, k). 
Hence, Lemma l2^ imr)lies that T 2 is a coordinatization. □ 

3. PFD-algorithm for Directed Hypergraphs 

3.1. Workflow 

We give here a summary of the workflow of the algorithm to compute the prime factor of connected directed 
hypergraphs. The top-level control structure is summarized in Algorithm[T]PFD_of _Di-Hypergraphs in which the 
subroutines Preprocessing (Alg.|2]i and Combine (Alg.|3]l are used. 


Algorithm 1 PFD of Di-Hypergraphs 
1; INPUT: A hypergraph H = (V, E); 

2: Ft <— Preprocessing!//); [>Now, H = (Vicx, £’/«) and vertex coordinates are known. 

3: Let Gaux = ({L ■ • ■, ^1,0) [>Here, k is the number of prime factors of [H} 2 . ; 

4: for each e 6 /?iex do 

5: Let / € |1,..., be the (unique) coordinate where distinct x, v 6 V(e) differ; 

6 : for ie/ = {l,...,A:j\ Lido 

7: if inc(e, i) ( /iiex then 

8: add edge ij to Gaux; 

9: end if 

10: end for 

11: end for 

12: Combine!//= (Flex, Flex), Gaux); 

13: OUTPUT: PFD ffj//, of H; 


As input of PFD_of _Di-Hypergraphs a connected hypergraph H - (V,E) is expected. First of all, subroutine 
Preprocessing is called. Here, the PFD of [H ]2 = and the respective coordinatization T of [H ]2 is computed 

by application of the algorithm of Imrich and Peterin 12111 . Then the vertices, the vertices within the arcs and the arcs 
are ordered in lexicographic order. This helps to achieve the desired time-complexity in later steps. 

By definition, coordinatization T of [H ]2 is a pre-coordinatization of H. By construction of T and Lemma lZ3l the 
pre-coordinatization T is at least as fine as the coordinatization of H w.r.t. its PFD. By Remarkl^it suffices to find a 
suitable partition of / = {1,..., k} to derive the prime factors of H. To this end, we initialize in Line|3]of Algorithm[T] 
the auxiliary graph Gaux, where each vertex i represents an element of 1. The edge set is left empty. We might later add 
edges in order trace back which equivalence classes of Ry have to be combined, i.e., all vertices within one connected 
components of Gaux will then be in one class of the respective partition of 1. 
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Algorithm 2 Preprocessing 
1: INPUT: A connected hypergraph H = (V, E); 

2: Compute D-FFD of [H]^ = and vertex-coordinates with the Imrich-Peterin-Algorithm iH; 

3; Compute the list of lexicographic ordered vertices (w.r.t. their coordinates); 

4: for each e e E do 

5: Reorder f(e) and h(e) w.r.t. the lexicographic order of the vertices; 

6: end for 

7: Compute the list E]sx of lexicographic ordered arcs; w.r.t. to the lexicographic ordered sets t(e) and then /;(«); 
8: return (Vicx, Uiex) with respective vertex-coordinates; 


Algorithm 3 Combine 

1: INPUT: A hypergraph H = (V, E) with pre-coordinates w.r.t the PFD of [//] 2 , a graph Gauxi 
2: Compute connected components of Gauxi 

3: for each e e £iex do 

4: Let j e |1,... ,fc) be the (unique) coordinate where distinct x,y e Vie) differ; 

5: Let /; be the connected component containing vertex j\ 

6: Assign color / to hyperedge e; 

7: end for 

8: compute [H }2 where edges ij obtains the unique color of e where i, j e Ffe); 

9: compute coordinates of all vertices in V in [El)^’, [>This pre-coordinatization is the PFD-coordinatization of El 

10: return D-PFD of H-, 


We continue to check in the for-loop in LineHIfTTlof Algorithm[T] if for each arc e e isiex that is contained in some 
y'-layer its “copies” are also contained in “incremental-neighboring” y-layers, i.e., we check if inc(e, i) e E\^x- If this 
is not the case, then we add the edge i j to Gaux- Finally, we use the information of the connected components I\,.. .,Ii 
of Gaux that partition the set I in order to determine the prime factors of the given hypergraph H. To this end, the 
subroutine Combine is called and an edge-colored 2-section is computed. That is, each edge e e £’([//] 2 ) that is 
contained in the copy of factor Gi with i € G obtains color s. In other words, all prime factors Gj, j £ G of [H ]2 are 
combined to a single factor of [H ]2 and the edges in the respective G-layers obtain color s. W.r.t. this coloring it is 
possible to efficiently determine new vertex coordinates in [H ]2 which is then a factorization coordinatization of H. 
We will show, that this leads to a “finest” coordinatization of H and hence, to the prime factors Hi,Hr of H. 

For illustrative examples see Figures [T] and |2l 

3.2. Correctness 

We are now in the position to prove the correctness of the algorithm PFD_of _Di-Hypergraphs, summarized in 
the following theorem. 

Theorem 3.1. Algorithm\T\is sound and complete. 

Proof. Given a hypergraph H - (V, E). We start with a preprocessing and call in Algorithm[T]the Algorithm|2l Here, 
the PFD of [//]2 = the respective coordinatization Y of [//]2 is computed. This coordinatization T is by 

definition a pre- coordinatization of H. Since [//]2 is an undirected graph, it is allowed to apply the algorithm of Imrich 
and Peterin 12111 . Finally, the vertices, the vertices within the arcs and the arcs are ordered in lexicographic order. The 
latter task is not important for the correctness of the algorithm, but for the time-complexity that we will consider later 
on. 

We are now in Line[3]of Algorithm[T] By construction of T and Lemma l231 the pre-coordinatization Y is at least as 
fine as the coordinatization of H w.r.t. its PFD. By Lemma l2.10l it suffices to find a suitable partition of / = {1,..., k}. 
To this end, we initialize in the auxiliary graph Gaux where each vertex i represents an element of /. The edge set is left 
empty. We might later add edges in order trace back which equivalence classes have to be combined, i.e., all vertices 
within one connected components of Gaux will then be in one class of the respective partition of /. 

Now consider the for-loop in LinelrlfTTI For each e e Fiex we check in which coordinates the vertices in Vie) 
differ. Since Y is a pre-coordinatization and by Lemma l2Gl this is exactly one coordinate for each hyperarc. Let 
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Figure 1: The prime hypergraph H (left-hand side) with arc set E{H) - 

{ei = ({11,12,13}, {12,13,14)),: ez = ({21,22,23), {22,23,24)), 63 = ({14), {12}), 64 

({ll},{ 21 )),e 5 = ({12},{22)),e6 = ({13,23),{13,23)),e7 = ({14},{24)),e8 = ({13),{23})} admits a non-trivial prime 
factorization of its 2-section (right-hand side) into fLjDAiz. The vertices of H are labeled w.r.t. its pre- 
coordinatization given by the coordinatization of [//Jz- For the first arc 64 of the lexicographic ordered arc set 
Flex = {e 4 , ei, es, eg, eg, 63 , 67 , ez) we can observe that the increment inc(e 4 ,2) = (({12}, {22)) = eg e E\^^. Analo¬ 
gously, inc(ei, 1) = ez, inc(e5,2) = eg, and inc(eg,2) = ez are all contained in Fiex- However, when we arrive at the 
hyperedge eg we obtain that inc(eg, 2) = ({14,24), {14,24}) is not an arc of E\^^. Hence, applying Algorithm[T] would 
lead to an edge 12 in Gaux, resulting in a connected auxiliary graph and, H would be determined as prime. A second 
example for an arc e e Fiex with inc(e, i) i Fiex is the edge e - e^. 



Figure 2: The non-prime hypergraph El - //iD/Zz is the product of the directed prime hypergraph H\ in Fig. [T] 
and an undirected hypergraph H 2 with two vertices and one hyperedge. The vertices of H are labeled w.r.t. its pre- 
coordinatization given by the coordinatization of [//]z. The auxiliary graph Gaux is initialized as the graph with 
three vertices and empty edge-set in Algorithm[T] While the increment inc(e, 2) = {(123,223), (123,223)} of the arc 
e - {(113,213), (113,213)} is still contained in £iex, the increment inc(e, 3) = {(124,224), (124,224)} is not. Hence, 
the edge 13 is added to Gaux- Since the increments of all hyperedges of the form {(/I j, i2j), (/I j, i2j)] are contained 
in Flex, no further edges will be added to Gaux- Hence, the sub-hypergraph induced by vertices with identical 2nd 
coordinate, i.e., the {1,3}-layers, constitute the copies the prime factor Hi, while the sub-hypergraphs induced by 
vertices with identical 1st and 3rd coordinate, i.e, the 2-layers, are copies of the prime factor //z. 



















e £ iiiex be a chosen arc and assume that all vertices in V(e) differ in the j-th coordinate. Now, it is checked if for arc 
e its “copies” are contained in each GJ-layer where w £ inc(V(e), i). If for some arc e £ we observe that there 
is no hyperedge inc(e, i) - (inc(f(e), /), mc(h{e), i) £ Giex then there is no “copy” of e in some y'-pre-layer through w 
with w £ inc{V{e), i). In this case we add the edge ij to Gaux if not already set. The latter tasks are repeated for all 
hyperarcs e £ £iex. 

Finally, in Line [12] the Algorithm |3] is called. The task of this subroutine is to combine the pre-coordinates and 
thus, the pre-layers in order to determine the layers of the final prime-factors of H. Let /i,..., L be the connected 
components of Gaux- Clearly, 11 = {/i,..., L) is a partition of I. Let each Ij having Ij elements. Lemma 12.101 and 
Remark|2]imply that T: V(H) —> XJ_[{1, .. .Jr] is a pre-coordinatization of H. It remains to show that 

(1.) T is a coordinatization and 

(2.) T is at least as fine as the coordinatization given by the PFD of H. 

Claim (L): By construction, all e £ Giex where the vertices differ in the /-th coordinate w.r.t. T are now contained 
in some L-layer where / £ L £ 11. Moreover, for all e in some L-layer the increments inc(e, j) with j + i and j i L 
must be contained in H, as otherwise we would have added the edge ij to Gaux and hence, j £ /j. As the latter is true 
for all L-layers contained in 11 we can apply Lemma l2.11l and conclude that T is a coordinatization of H. 

Claim (2.): Given the pre-coordinatization T of H. By construction of T and Lemma 1231 T is at least as fine as 
the coordinatization of H w.r.t. its PFD. Thus, there is a partition 11' = {/[,..., /'} of I w.r.t. the PFD of H. It remains 
to show that if there are two indices i, y £ L £ 11, then /, j are also contained in the same class of 11'. If i, j e C eTl 
then they are in same connected component Cj of Gaux- Hence, it suffices to consider pairs i, j e Is e Yl that are 
connected by an edge. Assume, for contradiction that i, j are in different classes of H'. W.l.o.g. let i £ /[ and j £ 7^. 
Moreover, let 77 = 77/j Hence, for all I £ I[ and thus, in particular for I - lit holds that for all arcs e in some 

fj-layer there is an arc (inc(f(e), /), mc{h{e), i) £ Tiiex. The same holds with the role of i and j switched. However, in 
this case we would not add the edge ij to Gaux, a contradiction. 

To finish the PFD-computation we have to compute T. To this end, we compute the 2-section [77] 2 with edges xy 
colored with color j whenever x and y are contained in some edge e that is contained in some y'-layer of 77. Lemma 
12.31 implies that T is also a coordinatization of [77 ]2 and hence, all edges with same color j in [77 ]2 are contained 
in the same equivalence class of 7 ?^([ 77 ] 2 ). In other words, (e,/) £ 7?^(77) if and only if (xy, uv) £ 7 ?^([ 77 ] 2 ) for all 
distinct x,y e V(e) and distinct u, v £ V(f). By construction, 7 ?;y([ 77 ] 2 ) is a a product relation of [ 77 ] 2 , and thus we 
can apply again a method proposed the by Imrich and Peterin (cf. Theorem 5.1. in iEUb . in order to obtain the desired 
coordinates and hence, T. □ 

3.3. Time Complexity 

In order to prove the time-complexity results, we first give the following lemma. 

Lemma 3.2. Let H be a hypergraph, let [77]2 = DjLjTf, be a factorization of its 2-section into k factors, and let m and 
n be the number of arcs and vertices ofH, respectively. Then for any I £ No it holds that log m — 0{n). 

Proof. If m, and u, are the numbers of arcs and vertices of the factors then m, < (2"')^, since we have to consider tail 
and head independently. Let N be the maximum number of vertices of a factor. Then we have 


k 


k 


k 


k 


k 



Taking logarithms on both sides of the inequality gives 


k 



for some suitable constant c. 
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On the other hand by bounding the size of every factor except the biggest one from below by 2 we get n - 
Y[]=i ni > N ■ 2^“^ Clearly, N ■ 2*^“' > c' ■ N ■ k' for some suitable constant c' depending on 1. Together with the 
estimate for log m this proves the lemma. □ 

The next two lemmas are concerned with the time-complexity of the subroutines Preprocessing and Combine. 

Lemma 3.3. Let H — (V, E) be a connected hypergraph with |y| = n, lEI = m and rank r. Then Algorithm\2\performs 
in O(r^mn) time. If we assume that H has bounded rank, then Algorithm^has time-complexity 0(m log^(«)). 

Proof. In Line|2] the first task is the computation of the 2-section [//] 2 - To this end, we initialize an adjacency 
list A^[l],.. .,N{n'\ with empty entries, which can be done in 0{n) time. We add for each arc e 6 £ and each pair 
{/, j] e the vertex i to £^[ 7 ] and j to A^[/], if these vertices are not already contained in the respective adjacency 
lists. Hence, we must check whether i e £^[ 7 ] or not. To this end, assume that £^[ 7 ] is already ordered. Hence we 
need 0{log{n)) comparisons to verify if i e N[j]. If this is not the case, the vertex i is added to £^[ 7 ] on the respective 
position so that £^[ 7 ] stays sorted. Analogously, we add 7 to A^[/], whenever 7 is not contained in A^[/]. As for each 
arc e e E there are at most ( 2 ) = O(r^) pairs {/, 7 } and for each such pair we have 0(log(n)) comparisons we end in 
a time-complexity of 0{n + mr^ log(«)) to create the adjacency list A[l],..., N[n]. These lists serve than as input for 
the algorithm of Imrich and Peterin which computes the PFD of the 2-section in 0{\E{{H'\2)\ + n) time. Since [//]2 
is connected and thus, [//]2 has at least |y| - 1 edges, the PFD algorithm runs in fact in 0{\E{{H'\2\) - 0{mr^) time. 
Hence, the total time complexity of Line|2]of Algorithm|2]is 0{n + mr^ log(u) + mr^) — 0{mr^ log(u)). 

In what follows, let k be the number of factors of [ 7 / 2 ] and let each v 6 y be identified with its respective (pre-) 
coordinate vector (vi,..., vf) computed by the Imrich-Peterin-Algorithm. Note, k is bounded by log(n). 

In Line 121 the list V of vertices is reordered in lexicographic order w.r.t. the vertex coordinates, i.e., v < w if there 
is some i e {1,..., k} with Vj < wj for all 7 e {1,..., / - 1} and v, < w,. This task can be done in 0{n log(n)k) = 
Ofi log^(n)). Since mr > n we obtain that 0(n log^(n)) = 0(mrlog^(n)). This new ordered vertex list is called Piex. 

We are now concerned with the for-loop in Line|4] For each hyperedge e e £ we reorder the vertices of its head 
and tail w.r.t. to the order of the vertices in Fiex. Each hyperedge contains at most r vertices and hence, this task can 
be done in 0(r\og{r)) time. Therefore, the entire for-loop (Line|4]-0 takes 0{mr\og{r)) time. 

Finally, the arcs are reordered w.r.t. the lexicographic ordered sets t{e) and h{e). We say e < / if t{e) < t(f) 
or t{e) = t(f) and h{e) < h{f), whereby the tails, resp., heads are compared w.r.t. the lexicographic order of their 
vertices. To determine if t(e) < t(f) or h{e) < h{f) for some arcs e,f e £, the at most 2r pairs of vertices must 
be compared, whereby the comparison of each such pair can be done in 0 ( 1 ) time, since the vertices are already 
ordered in the tails and heads. The reordering of the arcs need than 0{m \og{m)) comparisons, where each comparison 
can be done in 0{r) time, by the preceding arguments. Hence, the creation of £iex takes 0(rm\og(m)) time. By 
Lemma 122] this is 0{rmn). Moreover, if we assume that the rank r is bounded, then m < < rn'' = 0{n''). 

Hence 0(\og(m)) - 0{\og(n'')) - 0{r\og{n)) = 0(\og{n)). In this case the time complexity for determining £iex is 
0{m\og(m)) = 0{m\og{n)). 

Taken together the latter arguments, we end in overall time complexity for Algorithm|2]of O(r^mn) and if the rank 
r is bounded with 0(m log^(n)). □ 

Lemma 3.4. Let H — (V, E) be a connected hypergraph with |y| = n, |£| = m and rank r. Moreover, assume that the 
graph Gawc has k vertices and m' edge with k < log(n). Then Algorithm\^preforms in 0(mnr^) time. If we assume 
that H has bounded rank, then Algorithm^has time-complexity 0(m log^(n)). 

Proof. Determining the connected components of Gaux in Line|2]can be done in 0{k + m') - 0(log(n) ■+■ \og~(n)) time 
by application of the classical breadth-first search. While doing this, we will in addition record in G(l) time for each 
vertex in which connected component it is contained. Let /i,..., 7/ be the connected components of Gaux- 

For each of the m arcs we have to find the indices where the vertices of the particular arc differs. To this end, it 
suffices to take any two vertices x and y of V{e) and to compare their k coordinates which takes 0{k) time. Let 7 be 
the coordinate where the two vertices differ. We need to check in which of the connected components f the vertex 
is contained in, which can be done in G(l) time, since we have already recorded for each vertex of Gaux. in which 
component it is contained in. Now, the color for each arc can be recorded in G(l) time. Hence, the for-loop (Line[21l7]l 
has overall-time complexity 0{mk) — G(mlog(«)). 
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To compute the 2-section in Line|8]with colored edges we initialize an extended adjacency list A^[l],..., A^[n] 
where whenever we add some i e A^[y] we also record the respective unique color of i j as a 2nd parameter. Recording 
this parameter can be done in (9(1) time, as for each arc e e £ it is known which color it has. Hence, we can argue 
analogously as in the proof of Lemma l3.3l and state that the 2-section with additionally colored edges can be computed 
in 0{n + mp- log(n)) time. 

Finally, the vertex-coordinates in [H ]2 can be computed in O(m') - 0(mP) time, see Theorem 5.1. in iH. 

Hence the overall-time complexity of Algorithm [2 is (9(log^(n) + mlog(n) + n + mP\og(n) + mP) - 0(n + 
m\og^(n)P). Since mr > n and log^(n) = 0(n), the latter can be expressed as 0(mnP). If we assume in addition that 
the rank r is bounded we get 0{n -i- m log^(n)r^) = 0{mr + m log^(n)r^) = 0(m log^(n)). □ 

We are now in the position to determine the time-complexity of algorithm PFD_of _Di-Hypergraphs. 

Theorem 3.5. Let H — (V,E) be a connected hypergraph with |y| = n, lisl = m and rank r. Then Algorithm\I\ 
computes the PFD ofH in 0(mnr^) time. If the rank r is bounded the time-complexity of Algorithm\I\is 0(m log^(n)). 

Proof We suppose both the vertices and the hyperarcs of H implemented as integers and E implemented as an m x 2 
array, where each entry E[e, /] contains the list of vertices in t{e) if i = 1 and h{e) if i = 2. In Line |2] we call 
Preprocessing(//) which takes O(r^mn) time and if r is bounded 0(m\og^(n)) time (Lemma f3.31 l. 

In what follows, let k < log(«) be the number of factors of [H 2 ] and assume that each v e F is identified with its 
respective (pre-)coordinate vector (vi,..., vp. 

In Line [2 the auxiliary graph is be initialized. In particular, we initialize Gaux as adjacency list, i.e., we create 
empty lists A[l],..., A[k] which can be done in 0{k) time. 

We are now concerned with the for-loop in Line |4] - [TT] For each of the m arcs we have to find the indices 
where the vertices of the particular arc differs. To this end, any two vertices x and y of V{e) are chosen and their k 
coordinates are compared, which takes 0{k) time. The nested for-loop (Line l6] - fTOli is executed for all coordinates 
i where the vertices of arc e are identical and it is checked whether inc(e, /) = (inc(f(e), /), inc(/z(e), /) is contained in 
Giex or not. The increment (inc(f(e), (), inc(/i(e), i) can be computed in 0{r) time. Note, the vertices within inc(f(e), i) 
and inc(/z(e), i) are still lexicographically ordered as only vertex-coordinates are incremented that have been identical 
for the vertices within the arc and thus,their /-th positions are all still equal after the computation of inc(e, i). We now 
check whether (inc(f(e), i), inc(/i(e), i) i E\^^. Since Giex is already ordered, binary search finds the corresponding arc 
using at most 0(\og{m)) comparisons of arcs and since head and tail of each arc are in lexicographic order comparing 
two arcs takes 0{r) time. Therefore, the if-condition in Line Q takes 0{r -H rlog(m)) = G(rlog(OT)) time. In case, 
(inc(f(e), /), inc(/!(e),;) i Giex we have to add a respective edge ij in Gaux. if not already set. Hence, to check whether 
ij exists in Gaux. we need to validate if i 6 ^[J]. To this end, assume that is already ordered. Hence we need 
G(log(k)) comparisons to verify if i e ^[J]. If this is not the case i is added to on the respective position so that 
A[y] stays sorted. Similarly, j is added to A[(] whenever i i N[j]. Hence the nested for-loop in (Line l6] - fTOli has 
time complexity 0{k{r log(m) + log(k))) = 0{kr log(m)), since the number of arcs m is at least as big as the number 
(non-trivial) factors k. Take together the latter arguments, the entire for-loop in Line |4] - [TT] has time-complexity 
0{m{k + krlog(m))) = 0{mkrlog{m)). By Lemma l3T2l this is O(mnr). Moreover, if we assume that the rank r is 
bounded, then m <n^ and hence, G(log(m)) = G(log(«'')) = 0{r\og{n)) - G(log(u)) In this case, the time complexity 
of Linel4l-[TTIis 0{mk\og{m)) - 0{mk\og{n)) - 0{m\o^{n)). 

Finally, we use Algorithm |2 which performs in OfnnrP time and if the rank r is bounded it has time-complexity 
G(OTlog^(«)) (Lemma [3.4b . 

To summarize, each step of Algorithm [1] can be performed in 0{mnP) time and if the rank r is bounded the 
time-complexity is 0{m\o^{n)). □ 

Acknowledgment 

We thank the organizers of the 8th Slovenian Conference on Graph Theory (2015) in Kranjska Gora, where the 
authors participated, met and basically drafted the main ideas of this paper, while drinking a cold and tasty red Union, 
or was it a green Lasko? 

We also thank Wilfried Imrich and Iztok Peterin for helpful comments regarding the time-complexity of our 
algorithm. 


11 


References 


[1] Aurenhammer, F., Hagauer, J., Imrich, W., 1992. Cartesian graph factorization at logarithmic cost per edge. Comput. Complexity 2, 331-349. 

[2] Berge, C., 1989. Hypergraphs: Combinatorics of finite sets, volume 45. North-Holland, Amsterdam. 

[3] Black, T., 2015. Monotone properties of k-uniform hypergraphs are weakly evasive, in: Proceedings of the 2015 Conference on Innovations 

in Theoretical Computer Science, ACM, New York, NY, USA. pp. 383-391. 

[4] Bresar, B., 2004. On subgraphs of Cartesian product graphs and S-primeness. Discr. Math. 282, 43-52. 

[5] Bretto, A., 2006. Hypergraphs and the Helly property. Ars Comb. 78, 23-32. 

[6] Bretto, A., 2013. Applications of hypergraph theory: A brief overview, in: Hypergraph Theory. Springer International Publishing. Mathe¬ 
matical Engineering, pp. 111-116. 

[7] Bretto, A., Silvestre, Y, Vallee, T., 2009. Cartesian product of hypergraphs: properties and algorithms, in: 4th Athens Colloquium on 
Algorithms and Complexity (ACAC 2009), pp. 22-28. 

[8] Bretto, A., Silvestre, Y, Valle, T., 2013. Factorization of products of hypergraphs: Structure and algorithms. Theoretical Computer Science 
475,47-58. 

[9] Crespelle, C., Thierry, E., Lambert, T., 2013. A linear-time algorithm for computing the prime decomposition of a directed graph with regard 
to the Caitesian product, in: Du, D.Z., Zhang, G. (Eds.), Computing and Combinatorics. Springer Berlin Heidelberg, volume 7936 of Lecture 
Notes in Computer Science, pp. 469^80. 

[10] Dorfler, W., 1979. Multiple Covers of Hypergraphs. Annals of the New York Academy of Sciences 319, 169-176. 

[11] Feder, T., 1992. Product graph representations. J. Graph Theory 16, 467-488. 

[12] Feigenbaum, J., 1986. Directed Cartesian-product graphs have unique factorizations that can be computed in polynomial time. Discrete Appl. 
Math. 15, 105-110. 

[13] Feigenbaum, J., Hershberger, J., Schaffer, A., 1985. A polynomial time algorithm for finding the prime factors of Cartesian-product graphs. 
Discrete Appl. Math. 12, 123-138. 

[14] Hammack, R., Imrich, W., Klavzar, S., 2011. Handbook of Product Graphs. Discrete Mathematics and its Applications. 2nd ed., CRC Press. 

[15] Hellmuth, M., 2013. On the complexity of recognizing S-composite and S-prime graphs. Discrete Applied Mathematics 161, 1006 - 1013. 

[16] Hellmuth, M., Noll, M., Ostermeier, L., 2014. Strong products of hypergraphs: Unique prime factorization theorems and algorithms. Discrete 
Applied Mathematics 171, 60-71. 

[17] Hellmuth, M., Ostermeier, L., Stadler, R, 2012a. Diagonalized Cartesian products of S-prime graphs are S-prime. Discrete Math. 312, 74 - 
80. Algebraic Graph Theory - A Volume Dedicated to Gei1 Sabidussi on the Occasion of His 80th Birthday. 

[18] Hellmuth, M., Ostermeier, L., Stadler, P.F., 20l2b. A survey on hypergraph products. Math. Comput. Sci. 6, 1-32. 

[19] Imrich, W., 1967. Kartesisches Produkt von Mengensystemen und Graphen. Studia Sci. Math. Hungar. 2, 285 - 290. 

[20] Imrich, W., 1971. Uber das schwache Kartesische Produkt von Graphen. Journal of Combinatorial Theory 11, 1-16. 

[21] Imrich, W., Peterin, I., 2007. Recognizing Cartesian products in linear time. Discrete Math. 307, 472-483. 

[22] Kaveh, A., Alinejad, B., 2012. Hypergraph products for structural mechanics, in: Topping, B. (Ed.), Proceedings of the Eleventh International 
Conference on Computational Structures Technology, Civil-Comp Press, Stirlingshire, UK. doi 10.4203/ccp. 99.266 paper 266. 

[23] Kaveh, A., Alinejad, B., 2015. Hypergraph products for structural mechanics. Advances in Engineering Software 80, 72 - 81. Civil-Comp. 

[24] Klavzai', S., Lipovec, A., Petkovsek, M., 2002. On subgraphs of Cartesian product graphs. Discr. Math. 244, 223-230. 

[25] Ostermeier, L., Hellmuth, M., Stadler, PF., 2012. The Cartesian product of hypergraphs. Journal of Graph Theory 70, 180-196. 

[26] Sabidussi, G., 1960. Graph Multiplication. Mathematische Zeitschrift 72, 446^57. 

[27] Shapiro, H., 1953. The embedding of graphs in cubes and the design of sequential relay circuits. Bell Telephone Laboratories Memorandum 
Unpublished. 

[28] Sonntag, M., 1989. Hamiltonian properties of the Cartesian sum of hypergraphs. J. Inf. Process. Cybem. 25, 87-100. 

[29] Szamkolowisz, L., 1962. Remarks on the Cartesian product of two graphs. Colloq. Math. 9, 43—47. 

[30] Vizing, V.G., 1963. The Cartesian product of graphs. Vycisl. Sistemy No. 9, 30-43. 

[31] Winkler, R, 1987. Factoring a graph in polynomial time. European J. Combin. 8, 209-212. 

[32] Zhu, X., 1992. On the chromatic number of the product of hypergraphs. Ars Comb. 34, 25-31. 


12 


